Big O notation

In mathematics, big O notation is used to describe the limiting behavior of a function when the argument tends towards a particular value or infinity, usually in terms of simpler functions. It is a member of a larger family of notations that is called Landau notation, Bachmann-Landau notation, or asymptotic notation. In computer science, big O notation is used to classify algorithms by how they respond (i.e., in their processing time or working space requirements) to changes in input size.

Big O notation characterizes functions according to their growth rates: different functions with the same growth rate may be represented using the same O notation. A description of a function in terms of big O notation usually only provides an upper bound on the growth rate of the function. Associated with big O notation are several related notations, using the symbols o, Ω, ω, and Θ, to describe other kinds of bounds on asymptotic growth rates.

Big O notation is also used in many other fields to provide similar estimates.

Contents

Formal definition

Let f(x) and g(x) be two functions defined on some subset of the real numbers. One writes

f(x)=O(g(x))\mbox{ as }x\to\infty\,

if and only if there is a positive constant M such that for all sufficiently large values of x, f(x) is at most M multiplied by g(x) in absolute value. That is, f(x) = O(g(x)) if and only if there exists a positive real number M and a real number x0 such that

|f(x)| \le \; M |g(x)|\mbox{ for all }x>x_0.

In many contexts, the assumption that we are interested in the growth rate as the variable x goes to infinity is left unstated, and one writes more simply that f(x) = O(g(x)). The notation can also be used to describe the behavior of f near some real number a (often, a = 0): we say

f(x)=O(g(x))\mbox{ as }x\to a\,

if and only if there exist positive numbers δ and M such that

|f(x)| \le \; M |g(x)|\mbox{ for }|x - a| < \delta.

If g(x) is non-zero for values of x sufficiently close to a, both of these definitions can be unified using the limit superior:

f(x)=O(g(x))\mbox{ as }x \to a\,

if and only if

\limsup_{x\to a} \left|\frac{f(x)}{g(x)}\right| < \infty.

Example

In typical usage, the formal definition of O notation is not used directly; rather, the O notation for a function f(x) is derived by the following simplification rules:

For example, let f(x) = 6x^4 - 2x^3 %2B5, and suppose we wish to simplify this function, using O notation, to describe its growth rate as x approaches infinity. This function is the sum of three terms: 6x4, −2x3, and 5. Of these three terms, the one with the highest growth rate is the one with the largest exponent as a function of x, namely 6x4. Now one may apply the second rule: 6x4 is a product of 6 and x4 in which the first factor does not depend on x. Omitting this factor results in the simplified form x4. Thus, we say that f(x) is a big-oh of (x4) or mathematically we can write f(x) = O(x4). One may confirm this calculation using the formal definition: let f(x) = 6x4 − 2x3 + 5 and g(x) = x4. Applying the formal definition from above, the statement that f(x) = O(x4) is equivalent to its expansion,

|f(x)| \le \; M |g(x)|

for some suitable choice of x0 and M and for all x > x0. To prove this, let x0 = 1 and M = 13. Then, for all x > x0:

\begin{align}|6x^4 - 2x^3 %2B 5| &\le 6x^4 %2B |2x^3| %2B 5\\
                                      &\le 6x^4 %2B 2x^4 %2B 5x^4\\
                                      &\le 13x^4\\
                                      &\le 13|x^4|\end{align}

so

 |6x^4 - 2x^3 %2B 5| \le 13 \,|x^4 |.

Usage

Big O notation has two main areas of application. In mathematics, it is commonly used to describe how closely a finite series approximates a given function, especially in the case of a truncated Taylor series or asymptotic expansion. In computer science, it is useful in the analysis of algorithms. In both applications, the function g(x) appearing within the O(...) is typically chosen to be as simple as possible, omitting constant factors and lower order terms. There are two formally close, but noticeably different, usages of this notation: infinite asymptotics and infinitesimal asymptotics. This distinction is only in application and not in principle, however—the formal definition for the "big O" is the same for both cases, only with different limits for the function argument.

Infinite asymptotics

Big O notation is useful when analyzing algorithms for efficiency. For example, the time (or the number of steps) it takes to complete a problem of size n might be found to be T(n) = 4n2 − 2n + 2. As n grows large, the n2 term will come to dominate, so that all other terms can be neglected — for instance when n = 500, the term 4n2 is 1000 times as large as the 2n term. Ignoring the latter would have negligible effect on the expression's value for most purposes. Further, the coefficients become irrelevant if we compare to any other order of expression, such as an expression containing a term n3 or n4. Even if T(n) = 1,000,000n2, if U(n) = n3, the latter will always exceed the former once n grows larger than 1,000,000 (T(1,000,000) = 1,000,0003= U(1,000,000)). Additionally, the number of steps depends on the details of the machine model on which the algorithm runs, but different types of machines typically vary by only a constant factor in the number of steps needed to execute an algorithm. So the big O notation captures what remains: we write either

\ T(n)= O(n^2)

or

T(n)\in O(n^2)

and say that the algorithm has order of n2 time complexity. Note that "=" is not meant to express "is equal to" in its normal mathematical sense, but rather a more colloquial "is", so the second expression is technically accurate (see the "Equals sign" discussion below) while the first is a common abuse of notation.[1]

Infinitesimal asymptotics

Big O can also be used to describe the error term in an approximation to a mathematical function. The most significant terms are written explicitly, and then the least-significant terms are summarized in a single big O term. For example,

e^x=1%2Bx%2B\frac{x^2}{2}%2BO(x^3)\qquad\text{as } x\to 0

expresses the fact that the error, the difference \ e^x - (1 %2B x %2B x^2/2), is smaller in absolute value than some constant times |x^3| when x is close enough to 0.

Properties

If a function f(n) can be written as a finite sum of other functions, then the fastest growing one determines the order of f(n). For example

f(n) = 9 \log n %2B 5 (\log n)^3 %2B 3n^2 %2B 2n^3 = O(n^3)\,\!.

In particular, if a function may be bounded by a polynomial in n, then as n tends to infinity, one may disregard lower-order terms of the polynomial. O(nc) and O(cn) are very different. The latter grows much, much faster, no matter how big the constant c is (as long as it is greater than one). A function that grows faster than any power of n is called superpolynomial. One that grows more slowly than any exponential function of the form c^n is called subexponential. An algorithm can require time that is both superpolynomial and subexponential; examples of this include the fastest known algorithms for integer factorization. O(log n) is exactly the same as O(log(nc)). The logarithms differ only by a constant factor (since \log(n^c)=c \log n) and thus the big O notation ignores that. Similarly, logs with different constant bases are equivalent. Exponentials with different bases, on the other hand, are not of the same order. For example, 2^n and 3^n are not of the same order. Changing units may or may not affect the order of the resulting algorithm. Changing units is equivalent to multiplying the appropriate variable by a constant wherever it appears. For example, if an algorithm runs in the order of n2, replacing n by cn means the algorithm runs in the order of c^2n^2, and the big O notation ignores the constant c^2. This can be written as  c^2n^2 \in O(n^2) . If, however, an algorithm runs in the order of 2^n, replacing n with cn gives 2^{cn} = (2^c)^n. This is not equivalent to 2^n in general. Changing of variable may affect the order of the resulting algorithm. For example, if an algorithm's running time is O(n) when measured in terms of the number n of digits of an input number x, then its running time is O(log x) when measured as a function of the input number x itself, because n = Θ(log x).

Product

 f_1 \in O(g_1) \text{ and } f_2\in O(g_2)\, \Rightarrow f_1  f_2\in O(g_1  g_2)\,
f\cdot O(g) \subset O(f g)

Sum

 f_1 \in O(g_1) \text{ and }
  f_2\in O(g_2)\, \Rightarrow f_1 %2B f_2\in O(|g_1| %2B |g_2|)\,
This implies  f_1 \in O(g) \text{ and } f_2 \in O(g) \Rightarrow f_1%2Bf_2 \in O(g) , which means that O(g) is a convex cone.
If f and g are positive functions, f %2B O(g) \in O(f %2B g)

Multiplication by a constant

Let k be a constant. Then:
\ O(k g) = O(g) if k is nonzero.
f\in O(g) \Rightarrow kf\in O(g).

Multiple variables

Big O (and little o, and Ω...) can also be used with multiple variables. To define Big O formally for multiple variables, suppose f(\vec{x}) and g(\vec{x}) are two functions defined on some subset of \mathbb{R}^n. We say

f(\vec{x})\mbox{ is }O(g(\vec{x}))\mbox{ as }\vec{x}\to\infty

if and only if

\exists C\,\exists M>0\mbox{ such that } |f(\vec{x})| \le C |g(\vec{x})|\mbox{ for all }\vec{x} \mbox{ with } x_i>M \mbox{ for all }i.

For example, the statement

f(n,m) = n^2 %2B m^3 %2B O(n%2Bm) \mbox{ as } n,m\to\infty\,

asserts that there exist constants C and M such that

\forall n, m>M\colon |g(n,m)| \le C(n%2Bm),

where g(n,m) is defined by

f(n,m) = n^2 %2B m^3 %2B g(n,m).\,

Note that this definition allows all of the coordinates of \vec{x} to increase to infinity. In particular, the statement

f(n,m) = O(n^m) \mbox{ as } n,m\to\infty\,

(i.e., \exists C\,\exists M\,\forall n\,\forall m\dots) is quite different from

\forall m\colon f(n,m) = O(n^m) \mbox{ as } n\to\infty

(i.e., \forall m\,\exists C\,\exists M\,\forall n\dots).

Matters of notation

Equals sign

The statement "f(x) is O(g(x))" as defined above is usually written as f(x) = O(g(x)). Some consider this to be an abuse of notation, since the use of the equals sign could be misleading as it suggests a symmetry that this statement does not have. As de Bruijn says, O(x) = O(x2) is true but O(x2) = O(x) is not.[2] Knuth describes such statements as "one-way equalities", since if the sides could be reversed, "we could deduce ridiculous things like n = n2 from the identities n = O(n2) and n2 = O(n2)."[3] For these reasons, it would be more precise to use set notation and write f(x) ∈ O(g(x)), thinking of O(g(x)) as the class of all functions h(x) such that |h(x)| ≤ C|g(x)| for some constant C.[3] However, the use of the equals sign is customary. Knuth pointed out that "mathematicians customarily use the = sign as they use the word 'is' in English: Aristotle is a man, but a man isn't necessarily Aristotle."[4]

Other arithmetic operators

Big O notation can also be used in conjunction with other arithmetic operators in more complicated equations. For example, h(x) + O(f(x)) denotes the collection of functions having the growth of h(x) plus a part whose growth is limited to that of f(x). Thus,

g(x) = h(x) %2B O(f(x))\,

expresses the same as

g(x) - h(x) \in O(f(x))\,.

Example

Suppose an algorithm is being developed to operate on a set of n elements. Its developers are interested in finding a function T(n) that will express how long the algorithm will take to run (in some arbitrary measurement of time) in terms of the number of elements in the input set. The algorithm works by first calling a subroutine to sort the elements in the set and then perform its own operations. The sort has a known time complexity of O(n2), and after the subroutine runs the algorithm must take an additional 55n^3%2B2n%2B10 time before it terminates. Thus the overall time complexity of the algorithm can be expressed as

T(n)=O(n^2)%2B55n^3%2B2n%2B10.\

This can perhaps be most easily read by replacing O(n2) with "some function that grows asymptotically slower than n2 ". Again, this usage disregards some of the formal meaning of the "=" and "+" symbols, but it does allow one to use the big O notation as a kind of convenient placeholder.

Declaration of variables

Another feature of the notation, although less exceptional, is that function arguments may need to be inferred from the context when several variables are involved. The following two right-hand side big O notations have dramatically different meanings:

f(m) = O(m^n)\,,
g(n)\,\, = O(m^n)\,.

The first case states that f(m) exhibits polynomial growth, while the second, assuming m > 1, states that g(n) exhibits exponential growth. To avoid confusion, some authors use the notation

g(x) \in O(f(x))\,.

rather than the less explicit

g \in O(f)\,,

Multiple usages

In more complicated usage, O(...) can appear in different places in an equation, even several times on each side. For example, the following are true for n\to\infty

(n%2B1)^2 = n^2 %2B O(n)\
(n%2BO(n^{1/2}))(n %2B O(\log n))^2 = n^3 %2B O(n^{5/2})\
n^{O(1)} = O(e^n).\

The meaning of such statements is as follows: for any functions which satisfy each O(...) on the left side, there are some functions satisfying each O(...) on the right side, such that substituting all these functions into the equation makes the two sides equal. For example, the third equation above means: "For any function f(n)=O(1)\,, there is some function g(n)=O(e^n)\, such that n^{f(n)}=g(n)." In terms of the "set notation" above, the meaning is that the class of functions represented by the left side is a subset of the class of functions represented by the right side. In this use the "=" is a formal symbol that unlike the usual use of "=" is not a symmetric relation. Thus for example n^{O(1)} = O(e^n)\, does not imply the false statement O(e^n) = n^{O(1)}\, .

Orders of common functions

Here is a list of classes of functions that are commonly encountered when analyzing the running time of an algorithm. In each case, c is a constant and n increases without bound. The slower-growing functions are generally listed first.

Notation Name Example
O(1)\, constant Determining if a number is even or odd; using a constant-size lookup table or hash table
O(\log \log n)\, double logarithmic Finding an item using interpolation search in a sorted array of uniformly distributed values.
O(\log n)\, logarithmic Finding an item in a sorted array with a binary search or a balanced search tree as well as all operations in a Binomial heap.
O(n^c),\;0<c<1 fractional power Searching in a kd-tree
O(n)\, linear Finding an item in an unsorted list or a malformed tree (worst case) or in an unsorted array; Adding two n-bit integers by ripple carry.
O(n\log n)=O(\log n!)\, linearithmic, loglinear, or quasilinear Performing a Fast Fourier transform; heapsort, quicksort (best and average case), or merge sort
O(n^2)\, quadratic Multiplying two n-digit numbers by a simple algorithm; bubble sort (worst case or naive implementation), Shell sort, quicksort (worst case), selection sort or insertion sort
O(n^c),\;c>1 polynomial or algebraic Tree-adjoining grammar parsing; maximum matching for bipartite graphs
L_n[\alpha,c],\;0 < \alpha < 1=\,
e^{(c%2Bo(1)) (\ln n)^\alpha (\ln \ln n)^{1-\alpha}}
L-notation or sub-exponential Factoring a number using the quadratic sieve or number field sieve
O(c^n),\;c>1 exponential Finding the (exact) solution to the travelling salesman problem using dynamic programming; determining if two logical statements are equivalent using brute-force search
O(n!)\, factorial Solving the traveling salesman problem via brute-force search; generating all unrestricted permutations of a poset; finding the determinant with expansion by minors.

The statement f(n)=O(n!)\, is sometimes weakened to f(n)=O\left(n^n\right) to derive simpler formulas for asymptotic complexity. For any k>0 and c>0, O(n^c(\log n)^k) is a subset of O(n^{c%2B\varepsilon }) for any  \varepsilon >0, so may be considered as a polynomial with some bigger order.

Related asymptotic notations

Big O is the most commonly used asymptotic notation for comparing functions, although in many cases Big O may be replaced with Big Theta Θ for asymptotically tighter bounds. Here, we define some related notations in terms of Big O, progressing up to the family of Bachmann–Landau notations to which Big O notation belongs.

Little-o notation

The relation f(x) \in  o(g(x)) is read as "f(x) is little-o of g(x)". Intuitively, it means that g(x) grows much faster than f(x), or similarly, the growth of f(x) is nothing compared to that of g(x). It assumes that f and g are both functions of one variable. Formally, f(n) = o(g(n)) as n → ∞ means that for every positive constant ε there exists a constant N such that

|f(n)|\leq\epsilon|g(n)|\qquad\text{for all }n\geq N~.[3]

Note the difference between the earlier formal definition for the big-O notation, and the present definition of little-o: while the former has to be true for at least one constant M the latter must hold for every positive constant ε, however small.[1] In this way little-o notation makes a stronger statement than the corresponding big-O notation: every function that is little-o of g is also big-O of g, but not every function that is big-O g is also little-o of g (for instance g itself is not, unless it is identically zero near ∞).

If g(x) is nonzero, or at least becomes nonzero beyond a certain point, the relation f(x) = o(g(x)) is equivalent to

\lim_{x \to \infty}\frac{f(x)}{g(x)}=0.

For example,

Little-o notation is common in mathematics but rarer in computer science. In computer science the variable (and function value) is most often a natural number. In mathematics, the variable and function values are often real numbers. The following properties can be useful:

As with big O notation, the statement "f(x) is o(g(x))" is usually written as  f(x) = o(g(x)), which is a slight abuse of notation.

Family of Bachmann–Landau notations

Notation Name Intuition As  n \to \infty, eventually... Definition
f(n) \in O(g(n)) Big Omicron; Big O; Big Oh f is bounded above by g (up to constant factor) asymptotically |f(n)|  \leq  g(n)\cdot k for some k \exists k>0, \exists n_0, \; \forall n>n_0, \; |f(n)| \leq |g(n)\cdot k|
or
  \exists k>0, \exists n_0, \; \forall n>n_0, \; f(n) \leq g(n)\cdot k
f(n) \in \Omega(g(n)) (Note that, since the beginning of the 20th century, papers in number theory have been increasingly and widely using this notation in the weaker sense that f = o(g) is false) Big Omega f is bounded below by g (up to constant factor) asymptotically f(n)  \geq  g(n)\cdot k for some positive k \exists k>0, \exists n_0, \; \forall n>n_0, \; g(n)\cdot k \leq f(n)
f(n) \in \Theta(g(n)) Big Theta f is bounded both above and below by g asymptotically g(n)\cdot k_1 \leq f(n) \leq g(n)\cdot k_2 for some positive k1, k2 \exists k_1 >0 , \exists k_2 > 0, \exists n_0, \; \forall n>n_0,

g(n) \cdot k_1 \leq f(n) \leq g(n) \cdot k_2

f(n) \in o(g(n)) Small Omicron; Small O; Small Oh f is dominated by g asymptotically |f(n)| \le |g(n)|\cdot \varepsilon for every \varepsilon \forall \varepsilon>0, \; \exists n_0, \; \forall n>n_0, \; |f(n)| \le |g(n)\cdot \varepsilon|
f(n) \in \omega(g(n)) Small Omega f dominates g asymptotically f(n) \ge g(n)\cdot k for every k \forall k>0, \; \exists n_0, \; \forall n>n_0, \; g(n)\cdot k \le f(n)
f(n)\sim g(n)\! on the order of; "twiddles" f is equal to g asymptotically f(n)/g(n) \to 1 \forall \varepsilon > 0, \; \exists n_0, \; \forall n>n_0, \; \left|{f(n) \over g(n)}-1\right|<\varepsilon

Bachmann–Landau notation was designed around several mnemonics, as shown in the As  n \to \infty, eventually... column above and in the bullets below. To conceptually access these mnemonics, "omicron" can be read "o-micron" and "omega" can be read "o-mega". Also, the lower-case versus capitalization of the Greek letters in Bachmann–Landau notation is mnemonic.

Aside from Big O notation, the Big Theta Θ and Big Omega Ω notations are the two most often used in computer science; the Small Omega ω notation is rarely used in computer science.

Use in computer science

Informally, especially in computer science, the Big O notation often is permitted to be somewhat abused to describe an asymptotic tight bound where using Big Theta Θ notation might be more factually appropriate in a given context. For example, when considering a function T(n) = 73n^3 %2B 22n^2 %2B 58, all of the following are generally acceptable, but tightnesses of bound (i.e., bullets 2 and 3 below) are usually strongly preferred over laxness of bound (i.e., number 1 below).

  1. T(n) = O(n100), which is identical to T(n) ∈ O(n100)
  2. T(n) = O(n3), which is identical to T(n) ∈ O(n3)
  3. T(n) = Θ(n3), which is identical to T(n) ∈ Θ(n3).

The equivalent English statements are respectively:

  1. T(n) grows asymptotically no faster than n100
  2. T(n) grows asymptotically no faster than n3
  3. T(n) grows asymptotically as fast as n3.

So while all three statements are true, progressively more information is contained in each. In some fields, however, the Big O notation (number 2 in the lists above) would be used more commonly than the Big Theta notation (bullets number 3 in the lists above) because functions that grow more slowly are more desirable. For example, if T(n) represents the running time of a newly developed algorithm for input size n, the inventors and users of the algorithm might be more inclined to put an upper asymptotic bound on how long it will take to run without making an explicit statement about the lower asymptotic bound.

Extensions to the Bachmann–Landau notations

Another notation sometimes used in computer science is Õ (read soft-O): f(n) = Õ(g(n)) is shorthand for f(n) = O(g(n) logk g(n)) for some k. Essentially, it is Big O notation, ignoring logarithmic factors because the growth-rate effects of some other super-logarithmic function indicate a growth-rate explosion for large-sized input parameters that is more important to predicting bad run-time performance than the finer-point effects contributed by the logarithmic-growth factor(s). This notation is often used to obviate the "nitpicking" within growth-rates that are stated as too tightly bounded for the matters at hand (since logk n is always o(nε) for any constant k and any ε > 0). The L notation, defined as

L_n[\alpha,c]=O\left(e^{(c%2Bo(1))(\ln n)^\alpha(\ln\ln n)^{1-\alpha}}\right),

is convenient for functions that are between polynomial and exponential.

Generalizations and related usages

The generalization to functions taking values in any normed vector space is straightforward (replacing absolute values by norms), where f and g need not take their values in the same space. A generalization to functions g taking values in any topological group is also possible. The "limiting process" x→xo can also be generalized by introducing an arbitrary filter base, i.e. to directed nets f and g. The o notation can be used to define derivatives and differentiability in quite general spaces, and also (asymptotical) equivalence of functions,

 f\sim g \iff (f-g) \in o(g)

which is an equivalence relation and a more restrictive notion than the relationship "f is Θ(g)" from above. (It reduces to \lim f/g = 1 if f and g are positive real valued functions.) For example, 2x is Θ(x), but 2x − x is not o(x).

Graph theory

It is often useful to bound the running time of graph algorithms. Unlike most other computational problems, for a graph G = (V, E) there are two relevant parameters describing the size of the input: the number |V| of vertices in the graph and the number |E| of edges in the graph. Inside asymptotic notation (and only there), it is common to use the symbols V and E, when someone really means |V| and |E|. This convention simplifies asymptotic functions and make them easily readable. The symbols V and E are never used inside asymptotic notation with their literal meaning, since the number of vertices and edges must be non-negative, so this abuse of notation does not risk ambiguity. For example O(E %2B V \log V)\, means O((V,E) \mapsto |E| %2B |V|\cdot\log|V|)\, for a suitable metric of graphs. Another common convention—referring to the values |V| and |E| by the names n and m, respectively—sidesteps this ambiguity.

History

The notation was first introduced by number theorist Paul Bachmann in 1894, in the second volume of his book Analytische Zahlentheorie ("analytic number theory"), the first volume of which (not yet containing big O notation) was published in 1892.[5] The notation was popularized in the work of number theorist Edmund Landau; hence it is sometimes called a Landau symbol. It was popularized in computer science by Donald Knuth, who re-introduced the related Omega and Theta notations.[6] He also noted that the (then obscure) Omega notation had been introduced by Hardy and Littlewood[7] under a slightly different meaning, and proposed the current definition. Hardy's symbols were (in terms of the modern O notation)

 f\lesssim g \iff f \in O(g)   and    f\ll g \iff f\in o(g);

other similar symbols were sometimes used, such as \preceq and \prec\!\!\!\!\!\!\!\!\prec. The big-O, standing for "order of", was originally a capital omicron; today the identical-looking Latin capital letter O is used, but never the digit zero.

See also

References

  1. ^ a b Thomas H. Cormen et al., 2001, Introduction to Algorithms, Second Edition
  2. ^ N. G. de Bruijn (1958). Asymptotic Methods in Analysis. Amsterdam: North-Holland. pp. 5–7. ISBN 9780486642215. http://books.google.com/?id=_tnwmvHmVwMC&pg=PA5&vq=%22The+trouble+is%22. 
  3. ^ a b c Ronald Graham, Donald Knuth, and Oren Patashnik (1994). Concrete Mathematics (2 ed.). Reading, Massachusetts: Addison-Wesley. p. 446. ISBN 9780201558029. http://books.google.com/?id=pntQAAAAMAAJ&dq=editions:ISBN0201558025. 
  4. ^ Donald Knuth (June/July 1998). "Teach Calculus with Big O". Notices of the American Mathematical Society 45 (6): 687. http://www.ams.org/notices/199806/commentary.pdf.  (Unabridged version)
  5. ^ Nicholas J. Higham, Handbook of writing for the mathematical sciences, SIAM. ISBN 0-89871-420-6, p. 25
  6. ^ Donald Knuth. Big Omicron and big Omega and big Theta, ACM SIGACT News, Volume 8, Issue 2, 1976.
  7. ^ G. H. Hardy and J. E. Littlewood, Some problems of Diophantine approximation, Acta Mathematica 37 (1914), p. 225

Further reading

External links